Is the Unitarity of the quark-mixing-CKM-matrix violated in neutron /3-decay? 
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We report on a new measurement of neutron /3-decay asymmetry. From the result 
Ao = -0.1189(7), we derive the ratio of the axial vector to the vector coupling constant A = gA/gv 
= -1.2739(19). When included in the world average for the neutron lifetime r = 885.7(7)s, this gives 
the first element of the Cabibbo-Kobayashi-Maskawa (CKM) matrix V u d- With this value and the 
Particle Data Group values for V us and V u b, we find a deviation from the unitarity condition for the 
first row of the CKM matrix of A = 0.0083(28), which is 3.0 times the stated error. 
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As is well known, the quark eigenstates of the weak 
interaction do not correspond to the quark mass eigen- 
states. The weak eigenstates are related to the mass 
eigenstates in terms of a 3 x 3 unitary matrix V, the so 
called Cabibbo-Kobayashi-Maskawa (CKM) matrix. By 
convention, the u, c and t quarks are unmixed and all 
mixing is expressed via the CKM matrix V operating on 
d, s and b quarks. The values of individual matrix ele- 
ments are determined from weak decays of the relevant 
quarks. Unitarity requires that the sum of the squares of 
the matrix elements for each row and column be unity. 
So far precision tests of unitarity have only been possible 
for the first row of V , namely 
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In the Standard Model, the CKM matrix is unitary with 
A = 0. Usually, \V u d\ is derived from superallowed nu- 
clear /3-decay experiments to \V u d\ = 0.9740(5). This 
value includes nuclear structure effect corrections. Com- 
bined with kaon-, hyperon- and B-decays, this leads to 
A =0.0032(14), signaling a deviation from the unitarity 
condition by 2.3 a standard deviation 0]. However, some 
of the nuclear corrections are difficult to calculate, and 
therefore the Particle Data Group Q doubles the error 
in \V ud \- 

A violation of unitarity in the first row of the CKM 
matrix is a challenge to the three generation Standard 
Model. The data available so far do not preclude there 
being more than three generations; CKM matrix entries 
deduced from unitarity might be altered when the CKM 
matrix is expanded to accommodate more generations 
1^,1). A deviation A has been related to concepts be- 
yond the Standard Model, such as couplings to exotic 
fermions to the existence of an additional Z boson 

HQ or to the existence of right-handed currents in the 
weak interaction Q . Non-unitarity of the CKM matrix 
in models with an extended quark sector give rise to an 
induced neutron electric dipol moment that can be within 
reach of the next generation of experiments || . 

In this article, we derive IV^I, not from nuclear (3- 
decay, but from neutron decay data. In this way, the 
unitarity check of (1) is based solely on particle data, i.e. 
neutron /3-decay, K-decays, and B-decays, where theo- 
retical uncertainties are significantly smaller. So much 
progress has been made using highly polarized cold neu- 
tron beams with an improved detector setup that we are 
now capable of competing with nuclear /3-decays in ex- 
tracting a value for V u d, whilst avoiding the problems 
linked to nuclear structure. 

In the Standard Model, the Lagrangian of neutron de- 
cay is restricted to: 



C int = ^V ud -{pl^(l + \j 5 ) 



2m n 
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Gf is the Fermi decay constant; n, p, e and v are spinors 
describing neutron, proton, electron and neutrino; A is 



the ratio of the axial vector to the vector coupling con- 
stant gA/gv] and q is the momentum transfer between 
hadrons and leptons. The term (fi p — fx n )/2m p is the 
weak magnetism contribution, which is linked to fi p and 
the anomalous magnetic moments of proton and neu- 
tron. m v is the nucleon mass. In Eq. (2), we omit- 
ted the second-class induced scalar form factor and 
the induced pseudotensor form factor g 2 , because second 
class currents are excluded in the Standard Model. The 
first-class induced pseudoscalar form factor g 3 is negligi- 
ble in neutron decay and constrained by the Goldberger- 
Treiman relation. 

Since Gf is known from muon decay, in the Standard 
Model only two additional parameters are needed to de- 
scribe free neutron decay, namely A and V u d- In prin- 
ciple, the ratio A can be determined from QCD lattice 
gauge theory calculations, but the results of the best cal- 
culations vary by up to 30%. Today therefore all weak 
semileptonic particle cross-sections used in cosmology, as- 
trophysics and particle physics have to be calculated from 
neutron decay data. 

A neutron decays into a proton, an electron and an 
electron antineutrino. Observables are the neutron life- 
time r and spins tr e , ov, er p , and momenta p e , p v , p p of 
the electron, antineutrino and proton respectively. The 
electron spin, the proton spin and the antineutrino are 
not usually observed. The lifetime is given by 



= C|K d | 2 (l + 3A 2 )/^(l + A fl ), 
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where C = G F m 5 J{2ir 3 ) = 1.1613 ■ lO^s" 1 in h = c = 1 
units, f R =1.71482(15) is the phase space factor ( in- 
cluding the model independent radiative correction) ad- 
justed for the current value of the neutron-proton tran- 
sition energy. A R = 0.0240(8) [§(l|] is the model depen- 
dent radiative correction to the neutron decay rate, of 
which 0.0212 is straightforward electroweak-asymptotic 
QCD contribution, whereas the remaining 0.0028 de- 
pends on the strong interaction models. The neutron 
/3-decay rate and its relevant uncertainties at the 10~ 4 
level were reviewed recently |l2|,[l3| ■ 

The probability that an electron is emitted with an- 
gle f? with respect to the neutron spin polarization 



P = < a, > is 



W(d) = 1 + -PAcos(tf), 
c 
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where v /c is the electron velocity expressed in fractions of 
the speed of light. A is the /?- asymmetry coefficient which 
depends on A. On account of order 1% corrections for 
weak magnetism, gv~gA interference, and nucleon recoil, 
A has the form A = A (l+A fim (A 1 W + A 2 W + A 3 /W)) 
with electron total energy W = E e /m e c 2 + 1 (endpoint 
Wo). Aq is a function of A 



An 



, A(A + 1) 
' 1 + 3A 2 ' 
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where we have assumed that A is real. The coefficients 



Au m , A\, A 2l A 3 are from 



taking a different A 
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convention into consideration. In addition, a further 
small radiative correction of order 0.1% must be ap- 
plied. Other correlation coefficients (not measured in our 
experiment) are the antineutrino-electron correlation a, 
the antineutrino- asymmetry correlation B, and the time- 
reversal-violating triple correlation coefficient D. They 
also depend on A. Hence, various observables are acces- 
sible to experiment, so that the problem in extracting 
A and \V u d\ is overdetermined and, together with other 
data from particle and nuclear physics, many tests of the 
Standard Model become possible. Of course, a pertinent 
determination of the radiative corrections of Eq. (3) re- 
mains an important task. 

In the following section we report on our new mea- 
surement of the neutron /3-asymmetry coefficient A with 
the instrument PERKEO II, and on the consequences 
for | V u d | ■ Our first measurement fl6|| with this new spec- 
trometer gave a value for Aq which differed significantly 
from the combined previous data. In the light of this 
we decided to rcmcasure coefficient Aq with an improved 
setup. The new result confirms our earlier finding, with 
a reduced error but same result. PERKEO II was in- 
stalled at the PF1 cold neutron beam position at the 
High Flux Reactor at the Institut Laue-Langevin, Greno- 
ble. Cold neutrons are obtained from a 25 K deuterium 
cold moderator near the core of the 57 MW uranium re- 
actor. The neutrons are guided via a 60 m long neutron 
guide of cross-section 60 x 120 mm 2 to the experiment 
and are polarized by a supermirror polarizer of 30 x 45 
mm 2 cross section. The de Broglie wavelength spectrum 
of the cold neutron beam ranges from about 0.2 nm to 
1.3 nm. Above a wavelength of A > 1.3 nm, a very high 
degree of polarization is difficult to achieve. A long wave- 
length cut-off filter just in front of the supermirror po- 
larizer removes these undesired neutrons from the beam 
i pTf . The degree of neutron polarization was measured 
to be P = 98.9(3)% over the full cross-section of the 
beam. The polarization efficiency was monitored and it 
remained constant during the whole experiment. The 
neutron polarization is reversed periodically with a cur- 
rent sheet spin flipper, with measured spin flip efficiency 
of / = 99.7(1)%. 

The main component of the PERKEO II spectrometer 
is a superconducting 1 T magnet in a split pair configura- 
tion, with a coil diameter of about one meter. Neutrons 
pass through the spectrometer, whereas decay electrons 
are guided by the magnetic field to either one of two scin- 
tillation detectors with photomultiplicr readout. The de- 
tector's solid angle of acceptance is truly 2x27r above a 
threshold of 60 keV. Electron backscattering effects, se- 
rious sources of systematic error in /3-spectroscopy, are 
effectively suppressed. Technical details about the in- 
strument can be found in p8[ . 

The measured electron spectra Nj (E e ) and Nj;(E e ) 
in the two detectors (i=l,2) for neutron spin up and 
down, respectively, define the experimental asymmetry 
as a function of electron kinetic energy E e 



Nj(E e )-Nl(E e 
Nj(E e ) + Nl(E e 
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By using (4) and with < cos(tf) > = 1/2, A lexp (E) is 
directly related to the asymmetry parameter 

A exp (E e ) = A le „ p (E e ) - A 2 ^ p (E e ) = V -APf. (7) 

The main experimental errors are due to 
neutron spin polarization, background subtraction and 
detector response. To analyze the neutron spin 
polarization, a special setup of three additional spin 
flippers and two supermirror polarizers was used [19|. 
This gave an uncertainty of 0.3% in the measured beam 
polarization. As a very precise cross-check, in a sepa- 
rate setup, neutron polarization was measured again with 
three completely different methods: firstly, with our su- 
permirror setup, secondly, with a new method using an 
almost opaque 3 He spin filter |2(],|2l]], and thirdly, with 
a polarized proton filter p2| . The results of all three 
measurements agree to within 0.15%. 

A correction of 0.5% on A is due to background 
subtraction. One main feature of the PERKEO II spec- 
trometer is its high /3-decay count rate of about 270 s _1 
due to a large decay volume (80 x 80 x 270 mm 3 ) and a 
An detector. As a consequence, the signal- to-background 
rate in the range of interest (Fig. 1) was 7:1. Most of 
the background is environmental and was measured sep- 
arately and subtracted from the data. Extreme care is 
required to suppress any beam-related background, as 
discussed in [|16| . The beam divergence was limited to 
12 mrad by appropriate neutron baffles upstream, made 
from enriched 6 LiF plates. The beamstop was positioned 
6 m downstream of the decay volume. The /3-detectors 
were installed far off the beam at a transverse distance 
of 960 mm, and had no direct view to the polarizer, the 
baffles or the beam stop. Any beam-related gamma quan- 
tum had to undergo multiple directional changes before 
reaching the detector. Compared with our previous mea- 
surement with the same apparatus ||l6|| , this beam- related 
background was reduced by a factor of three to 0.3 s _1 
or to 1:200 of the signal rate. Thus, in the fit inter- 
val the size of the background correction is 0.5%. This 
beam-related background was determined using an ex- 
trapolation procedure described in Q and [ p^[ . The 
relative uncertainty of 50% is a conservative estimate of 
this background extrapolation method (see Table 1). 

The detector response function was determined with 
six conversion line sources on 10/ig/cm 2 carbon backings, 
which were remotely inserted into the spectrometer. The 
K, L, M and N conversion electrons and the correspond- 
ing Auger electrons are taken into account. Detector re- 
sponse is linear in energy within 1% leading to an uncer- 
tainty of 0.2% in A ||. 

The experimental function Ai (E e ) and a fit with one 
free parameter Ai^ v (the absolute scale of ^4o) is shown 
in Fig. |l|. The x 2 is 142.5 for detector one and \ 2 is 129 
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for detector two for 150 degree of freedom. The fit inter- 
val was chosen such that the signal-to-background ratio 
was at maximum. The parameter Ai cxp is directly related 
to the asymmetry parameter via Ai^ p = A^-P-f. From 
the experimental asymmetries we get Ay = -0.1174(7) 
and A2 exp = -0.1163(7) for detector 1 and detector 2 re- 
spectively. All corrections and uncertainties entering the 
determination of A are listed in Tab. 1. The corrections 
marked with an asterisk are already included in the fit. 
After correcting for the other small experimental system- 
atic effects listed in Tab. 1, we obtain Aq = -0.1189(8). 
This value is identical to our earlier result (li| of Aq — 
-0.1189(12), but with a smaller error. The combined re- 
sult is 



A a = -0.1189(7) and A = -1.2739(19). 



(8) 



With this value, and the world average value for r = 
885.7(7) s from §, we find from (3) that 



\V ud \ = 0.9713(13). 



(9) 



With j§ \V US \ = 0.2196(23) and the negligibly small \V ub \ 
= 0.0036(9), we obtain 



\V ud 



\Vus 



\V ub \ 2 = 1 - A = 0.9917(28). (10) 



This value differs from the Standard Model prediction by 
deviation A = 0.0083(28), or 3.0 times the stated error. 

Earlier experiments |22-|28| gave significantly lower 
values for |A|. However, in all these earlier experiments 
large corrections had to be made for neutron polarization, 
electron-magnetic mirror effects or background, which 
were all in the 15% to 30% range. In our experiment 
on the other hand, the total correction to the raw data 
is 2.0%, i.e. 10 times less than in earlier experiments. 
We therefore believe that our new experiment is more 
reliable than previous experiments. 

Averaging over our new result and previous neutron 
/3-decay results the Particle Data Group |23|] arrives at a 
new world average for |V^<j| from neutron /3-decay which 
leads to only a 2.2 a deviation from unitarity. 

The Particle Data Group obtains from superallowed 
0+ -> 0+ transitions a value \V ud \ = 0.9740(10). This 
value is compatible with our value (9) at the 90 % C.L.. 

An independent test of CKM unitarity comes from W 
physics at LEP 24 where W decay hadronic branching 
ratios can be used. Since decay into the top quark chan- 
nel is forbidden by energy conservation one would expect 
J2 \ Vij\ 2 to be 2 with a three generation unitary CKM 
matrix. The experimental result is 2.032(32), consistent 
with (10) but with considerably lower accuracy. 

In the frame of the present article we do not want to 
speculate on the origin of the deviation A. Nevertheless, 
we want to point out that it is unlikely that this devia- 
tion is due to induced form factors (as discussed above) 
or erroneous radiative corrections or to the other CKM 
elements V us and V u b- A non-zero second-class term gi 



(in contradiction to the Standart Model) is unlikely be- 
cause the present experimental limit on g2 < 0.2 |25} | 
would lead to a change in the neutron decay-asymmetry 
Aq by less than 0.15%. If the deviation A was due to the 
radiative correction A#, than the error on A# must be 
more than 9 times larger than the quoted error. Also it 
is unlikely that a nonzero A is due to an error in the de- 
termination of the high energy results on V us because the 
error in V us must be enlarged by factors of 8 to explain 
our value of A (the value of V u b is completely negligible 
in this context). 

In summary, \V u d\, the first element of the CKM ma- 
trix, has been derived from neutron decay experiments 
in such a way that a unitarity test of the CKM matrix 
can be performed based solely on particle physics data. 
With this value, we find a 3 a standard deviation from 
unitarity, which conflicts the prediction of the Standard 
Model of particle physics. 

This work was funded by the German Federal Min- 
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FIG. 1. Fit to the experimental asymmetry A exp for de- 
tector 1 and detector 2. The solid line shows the fit interval, 
whereas the dotted line shows an extrapolation to higher and 
lower energies. 
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Table 1: Experimental corrections and uncertainties 
entering the determination of A. 



Effect 


Correction 


Error 


Polarization analysis 






Polarization efficiency 


1.1% 


0.3% 


spinflip efficiency 


0.3% 


0.1% 


P)ata set 






Statistics 




0.45% 


Background 


0.5%* 


0.25% 


L/CLcLLUl 1 Co UUllftC 






Linearity 




0.2% 


Width and pedestal 




0.1% 


Drifts 




0.06% 


Edge effect 


-0.24%* 


0.1% 


Hemisphere integration 






Mirror effect 


0.09% 


0.02% 


backscattering 


0.2% 


0.17% 


Radiation corrections 


0.09% 


0.05% 


Sum 


2.04% 


0.68% 



* already included in the fit function of Fig. 1. 
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